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ABSTRACT 


Non -empirical  calculatione  of  the  electronic 
structure  of  the  ammonia  molecule  are  carried  out  uaing 
a  ''bond -orbital'*  and  a  "modified  electron-pair"  function. 
The  calculated  binding  energy  is  0.  348  a.u.  (76  o/o  of  the 
experimental  value)  and  0.381  a.u.  (83  o/o)  respectively. 
The  best  theoretical  value  previously  known  was  0.378  a.u. 
the  resiilt  of  an  self-consistent  LCAO  calculation  including 
limited  configuration  interaction. 

It  is  concluded  that  the  "modified  electron-petir" 
function  is  an  adequate  simple  approximation,  and  the 
"frozen  core"  approximation  employed  in  the  calculation 
is  satisfactory. 

The  dipole  moment  is  also  calculated,  and 
reasonable  agreement  with  the  experimental  value  is 
obtained. 
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l.  INTRODUCTION 

Hurleyt  Lennard-Jonet  and  Pople  and  Slater  pointed  out  the  merit 
of  using  a  set  of  orbitals  of  the  following  nature  in  the  calculation  of  the  elec* 
tronic  structure  of  molecules*  Two  orbitals  which  describe  one  and  the  same 
bond  are  not  orthogonal*  Orbitals  describing  different  bonds  are  orthogonal* 

Thus  there  is  no  difficulty  connected  with  non*orthogonality»  and  energy  expres¬ 
sions  can  be  given  in  a  fairly  simple  closed  fo'  m*  Such  orbitals  may  be  said 
to  have  pair  non-orthogonality* 

The  idea  was  further  developed  by  Parks  and  Parr  ^^*  However*  these 
authors  suggested  no  method  of  constructing  such  orbitals  with  pair  non- 
-orthogonality*  McWeeny  and  Ohno  proposed  one  way  of  doing  this*  They  first 
obtained  a  set  of  orthogonal  orbitals  and  then  made  the  orbitals  within  each  pair 
non -orthogonal  to  each  other*  They  carried  out  calculations  of  the  water 
molecid^e  using  these  orbitals*  In  one  of  their  calculations,  they  utilized  a 
^'modified  electron-pair'*  fimction  which  has  the  form  of  an  ordinary  electron- 
-pair-bond  function  but  is  constructed  from  orbitals  with  pair  non-orthogonality* 
In  this  method*  the  water  molecule  was  treated  effectively  as  a  four-electron 
problem,  and  the  wave  function  contained  only  two  parameters*  The  calculated 
total  energy  was  lower  than  that  given  by  the  elaborate  self-consistent  LCAO 
method  (the  SC  LCAO  method)  . 

In  the  modified  electron -pair  function*  McWeeny  and  Ohno  used  pre- 
-determined  orbitals  to  describe  the  inner-shell  and  the  lone  pair*  This  was 
called  the  "frozen  core"  approximation,  and  it  turned  out  to  be  quite  satisfactory* 

Being  encouraged  by  these  results  in  the  case  of  water*  we  apply  the 

modified  electron -pair  method  to  the  ammonia  molecule  in  this  paper*  We  adopt 

7) 

the  Hartree-Fock  atomic  orbitals  and  integral  values  calculated  by  Kaplan 
It  is  interesting  to  see  whether  the  differences  in  the  number  of  lone  pairs  (one 
for  two  for  H2O)  and  in  the  basic  orbitals  (Hartree-Fock  for  NH^,  Slater 

for  H2O)  affect  the  general  features  of  the  result* 


It  was  noted  by  Kaplan  that  changing  the  atomic  orbitals  from  Slater  to 
Hartree-Fock  orbitals  causes  appreciable  change  in  the  value  of  many 
molecular  integrals  ^)* 
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A  pertinent  remark  on  McWeeny  and  Ohno'e  paper  was  made  by  Merri* 
field  He  pointed  out  that  the  work  is  based  on  approximate  values  for  three - 
and  four-centre  integrals  and  on  single-term  Slater  orbitals.  He  seems  to 
feel  that  errors  caused  by  these  two  approximations  balance  one  another  to 
produce  fortuitous  results.  In  the  present  ammonia  calculation  we  use  Kaplan's 
integral  values  which  are  correct  within  +  10*^  a.u.  It  will  be  interesting  to 
see  whether  or  not  this  improvement  substantiates  Merrifield's  criticisms. 

We  carsry  out  an  additional  calculation  using  a  "bond -orbital"  function, 
as  the  bond -orbital  function  is  a  special  case  of  the  "modified  electron-pair" 
function. 


2.  ORBITALS  AND  INTEGRALS 

The  methods  of  obtaining  the  orthonormal  hybridised  orbitals  and  of 
transforming  the  molecular  integrals  to  integrals  of  these  new  basic  orbitals 
are  quite  similar  to  the  case  of  water  Therefore,  we  shall  only  briefly 
sketch  it  here. 

The  co-ordinate  system  to  be  used  is  shown  in  figure  1. 

We  start  with  the  set  of  orbitals  used  by  Kaplan: 

cl=(ny6Zx.(^  H„  Hx  ),  p.i) 

where,  n,  s,  z,  x,  and  y  are  the  Hartree-Fock  is,  28,  2pz,  2px  and  2py 

orbitals  of  the  nitrogen  atom  respectively.  The  orbitals  H  •  H  and  H  were 

y 

defined  as  linear  combinations  of  three  hydrogen  Is  orbitals  h^,  1^2$  and  h^  : 

H.  =  {3( 

H,  =  {  3(  I  -  Sk)4  \ 

It  is  convenient  to  go  back  to  the  orbitals  h^,  h2  and  h^  instead  of  H^, 
H^  and  Hy  at  this  stage: 
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L=('KxiBX^  i,,  -^1  -^3  )=  ,  <2.3) 

where  the  transformation  is  inverse  to  the  transformation  (2.2)* 

Since  the  nitrogen  is  electrons  have  much  lower  energy  than  the 
others*  we  shall  retain  the  nitrogen  Is  orbital  unchanged.  The  orbitals  h^> 
and  h^  are  made  orthogonal  to  the  nitrogen  is  orbital  by  the  Schmidt 
procedure.  This  yields  the  set* 

c=  ( 'H  -6  z  X  ^  fi"  -kl  ^  (2.4) 

We  treat  all  the  valence  orbitals*  28*  2p  of  N  and  Is  of  H  •  on  an 

9) 

equal  basis  and  apply  Ldwdin's  symmetrical  orthogonalization 

d  =  (o^^i  =  ^Sc  (2-6) 

where  Sq  overlap  matrix  for  the  orbitals  C  .  This  orthogonalization 

procedure  has  the  advantage  that  the  resulting  orbitals  resemble  the 

2 

original  orbitals  O  most  closely  in  the  sense  ^/|d.  -  c.  |  dv 
,  ,  .  10)  1  1  i 

is  a  minimum  '« 

The  final  step  is  to  get  s  *  z  >  x  ,  and  y  hybridized.  This  hybridiza* 
tion  is  carried  out  so  that  the  resulting  bond  hybrids  b^  ,  b2  and  b^  point 
along  the  bond  directions: 

e  =  (  -n  A  b,  bj  bj  I;  )  =  d  17.  (2.7) 

The  transformation  U  ha.  the  following  form, 

^  °  ° .  0  000 
0  -t  ->Jl  -Jl  -tV'VS  0  0  0 

0  si  I  -t2  -t/«J3  -t/Nfs  -t/sf5  006 

0  0  sji/sji  -l/fjh  0  0  0 

0  0  0  l/sjl  -  l/sli  0  0  0 

0  0  0  0  0  1  0  0 

0  0  0  0  0  0  1  0 

0  0  0  0  0  0  0  1 
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where  t  *  cot  0  ,  and  0  is  the  angle  between  a  bond  direction  N>H  and 
the  -Z  axis.  In  the  equilibrium  case,  ^HNH  is  106^47'.  and  cot  0  a 
>  0.404  865. 

The  final  orthonormal  hybridized  orbitals  are  related  to  Kaplan's 
basic  set  by 

e  =  (2.9) 

where  V  is  completely  determined  by  the  shape  of  the  molecule  and  the 
overlap  integrals  between  the  basic  atomic  orbitals# 

The  trainsformation  of  all  integrals  to  the  new  basis  G  can  be  carried 

4) 

out  in  the  same  way  as  in  the  case  of  water 

The  numerical  calculation  is  performed  for  the  equilibrium  case^  the 
N*H  distance  being  1. 9l6  a.u.  All  integrals  are  taken  from  Kaplan's  paper 
except  for  the  following 

(ss|sH^)  ^  -0#  1193  instead  of  -0.0864  , 

(ns|H^H^)  =  0.0001  instead  of  0.0014  , 

(sH^lzHj  =  0.0128  instead  of  0.0060  , 

(n|  r^|z  )  s  0.0070  instead  of  0*0076  • 


3.  THE  BOND-ORBITAL  APPROXIMATION 

The  simplest  electronic  wave  fxinction  is  obtained  in  the  following  way. 
Two  electrons  are  put  into  the  fixed  orbitals  n  and  i  in  order  to  represent 
the  inner-shell  and  the  lone  pair.  Each  bond  is  described  by  using  a  bond  orbital 
which  is  a  linear  combination  of  b^^  and  h.^  » 

Bi  =  (  +  t=  1,2.3.  (3.1) 

The  total  wave  function  is  then  a  single  Slater  determinant^ 

B,') 

=  ^  |'KCI)U(2)  i(3)  •  '  •  B»(8)B»(9)Bi,Cl0)|  ^ 


(3.2) 
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where  the  barred  orbitals  are  occupied  by  an  electron  with  spin  p  • 

Since  the  five  orbitals  in  (3«2)  are  orthogonal  to  each  otheri  the  total 
electronic  energy  is  expressed  by  a  simple  formula  which  contains  only  one 
variational  parameter  \  • 

2 

The  total  energy  reaches  a  minimum  when  X  a  0*42  ,  and  the  binding 

energy  is  0*348  a*u*  a  9.47  eV  which  is  96  o/o  of  the  binding  energy  calculated 

by  the  SC  LCAO  method.  The  experimental  binding  energy  is  0*459  a*u«  a 

2 

s  12*49  eV  *  The  effect  of  variation  of  X  on  the  energy  is  illustrated  in 
figure  2. 

2 

The  calculated  dipole  moment  for  X  a  0*42  is  1*84  Debye  while  that 
by  the  SC  LCAO  method  is  1*82  Debye  and  the  experimental  value  is  1*46  Debye. 

It  is  remarkable  that  this  simple  wave  function  produces  almost  as  good 

results  as  the  elaborate  SC  LCAO  function.  It  is  also  interesting  to  note  that 

2 

the  optimum  value  of  X  *  which  is  a  measure  of  the  bond  polarity,  is  very 
close  to  the  value  (0*43)  for  the  O-H  bond  in  water. 

Kaplan  constructed  his  equivalent  orbitals  from  his  SC  LCAO  molecular 
orbitals  in  the  following  way, 

It'  =  K 
i'=  . 

where  .  t|<Q  ,  .  •('jy  molecular  orbitals  in  the  order  of 

increasing  orbital  energy. 

We  may  compare  our  orbitals  n  ,  i  ,  B  corresponding  to  the  inner  shell, 
lone  pair  and  bond  orbitals  respectively  with  Kaplan's  equivalent  orbitals.  In 
terms  of  the  basic  set  ^  ,  they  are  expressed  as 

n  =  n  ^ 

i  =  0.020  n  -  0.700  s  +  0.781  z  -  0.210  , 

B,»  -0.025  n  -  0.292  s  -  0.231  z  +  0.530x  +  0,249  H  +  0.371  H 
1  O  ‘  .  x> 


The  choice  of  pure  as  the  lone*pair  orbital  is  rather  arbitrary.  In 

the  Cj^  case,  we  cannot  determine  equivalent  orbitals  by  symmetry 
consideration  only* 
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and 


n' 

i' 


S 


1.000  n  -  0.003  •  -  0.002  z  •  0.001  H 

o; 

0.026  n  -  0.442  8  4  0.896  z  -  0.258  H  . 

o ' 

-0.017  n  -  0.438  $  -  0.093  z  +  0.560x  +  0.157  H  +  0.397  H 

o 


(3.5) 


We  can  zee  a  marked  resemblance  between  these  two  sets  of  orbitals  although 
the  methods  of  construction  are  entirely  different. 


4.  THE  MODIFIED  ELECTRON -PAIR  APPROXIMATION 

As  for  the  inner-shell  and  the  lone  pair,  we  use  the  same  "frozen-core" 
approximation  as  in  the  previous  section.  The  description  of  the  three  bonds  is 
improved  by  using  the  following  two  functions  for  a  bond  instead  of  just  B.  , 

^ 

i  =  1,  2,  3  .  (4.1) 

We  use  the  electron-pair  function 

to  describe  a  bond  and  thus  the  total  wave  function  is  written  in  an  abbreviated 
form  as 

Jl'  [bt?,]  C  (4.3) 

The  orbitals  ET.  and  £!  are  not  orthogonal  any  more*  and  this  non- 
-orthogonality  within  a  bond  is  necessary  for  a  single  spin-paired  "structure" 
to  be  a  reasonably  good  approximation*  HoweveXi  the  non-orthogonalities  are 
limited  within  bonds,  and  we  can  still  express  the  total  electronic  energy  in  a 
fairly  simple  closed  form 

The  function  contains  two  variational  parameters  X  and  X'  »  and  for 
X  a  X'  the  present  approximation  reduces  to  the  bond  orbital  approximation. 

The  calculated  results  are  shown  in  Table  I*  The  dependence  of  the 


TABLE  I. 


Total  and  binding  energies  and  dipole  moment  in  various  approx¬ 
imations. 


Total  energy 
(a.u.) 

Binding  energy 
(a.u.) 

Dipole  moment 
(Debye) 

B<^d  orbital 
=  0.4Z) 

-56.250 

0.348 

1.84 

Modified  electron  pair 
(X2  =  0.20.  X  2  =  0.6^ 

-56.283 

0.361 

1.89 

SC  LCAO 

j 

-56.266 

0.364 

1.82 

SC  LCAO  +  Cl 

-56.280 

0.378 

- 

7) 

Experimental  ' 

1 

-56.  595 

' 

0.4S9 

1.46 

2  ’2 

energy  on  X  and  X  is  illustrated  in  figure  3* 

The  resulting  binding  energy  is  0.381  a.u.  s  10.37  eV  which  is  0.90  r-v 
better  than  the  bond  orbital  value  and  is  83  o/o  of  the  experimental  value. 

This  simple  treatment*  in  which  we  reduce  the  problem  effectively  to  that 

of  six -electrons,  gives  an  energy  0.46  eV  better  than  the  ten-electron  SC  LCAO 

2  *2 

calculation*  The  optimum  values  of  the  parameters  X  and  X  are  again  not 
so  different  from  those  for  water*  the  latter  being  0.  17  and  0.72. 

It  is  interesting  to  note  that  the  energy  value  calculated  with  the  modified 
electron-pair  function  is  slightly  better  than  that  obtained  by  the  limited  con¬ 
figuration  interaction  calculation  In  the  latter,  thirteen  configurations  were 
carefully  selected  from  the  total  ninety -eight  configurations. 

The  modified  electron-pair  method  seems  to  be  a  nice  compromise 
between  accuracy  and  ease  of  use. 
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5.  DISCUSSION 

We  have  used  orthogonalized  atomic  orbitals  as  our  basis.  The  aim  of 
this  was  not  only  to  circumvent  the  non  •orthogonality  difficulty  but  also  to 
emphasize  the  localization  of  bonds  and  lone  pairs.  Whether  the  present  basic 
orbitals  (2.7)  are  localized  or  not  is  difficult  to  say.  For  example*  the 
orbital  h^^  is  expressed  as 

-0,080  'l^  +  0.3^9  /6  +  O.ll'l  z  -0.351  X  4- 

+  1.345  ^,-0.123^^-0.123^3* 

This  looks  more  delocalized  than  the  original  hydrogen  orbital  h^  at  the  first 
sight.  However*  the  coefficients  of  the  other  atomic  orbitals  are  such  that  they 
reduce  the  electron  density  near  the  nitrogen  atom  and  hydrogen  atoms  2  and  3, 
so  that  the  orbital  h^  may  well  be  more  localized  than  h^  .  We  cannot 
compare  the  degree  of  localization  unless  we  introduce  a  certain  measure. 

Hall  suggested  the  use  of  self-Coulomb  integral  (  9P^|<P9)  for  this  measure 
The  values  for  and  h^  are 

( c  1  t'C) 

(  ^  ft.ll  t.ft,.  ) 

This  criterion*  therefore*  shows  that  h^^  is  the  one  more  localized. 

The  modified  electron-pair  function  is  satisfactory  in  that  it  gives  lower 
total  energy  than  either  the  SC  LCAO  method  or  the  limited  configuration 
interaction  method.  Thus  the  present  orthogonalization  procedure  to  keep  the 
nitrogen  Is  orbital  unperturbed  is  justified.  To  describe  the  inner  shell  and 
lone  pair  by  the  fixed  orbitals  n  and  i  respectively  -  this  may  be  called  the 
"frozen-core"  approximation  -  proves  also  satisfactory. 

The  wave  function  furnishes*  at  the  same  time*  a  clear-cut  picture  of 
bonding  and  lone -pair  electrons. 

In  the  electron-pair  function*  correlation  between  the  two  electrons  with 


=  0.  6899  a.u. 
=  0. 6250  a.u. 


(5.2) 


It  should  be  noted  that  the  sign  of  the  2s  orbital  of  nitrogen  is  opposite 
to  the  usual  one. 
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spin  Q  and  p  in  each  bond  is  taken  into  account  to  some  extent*  On  the  other 

hand,  this  correlation  is  not  at  all  taken  into  account  in  the  bond  orbital  func* 

Z  *2 

tion*  The  optimum  values  of  \  and  X  differ  considerably  from  that  in 
the  bond  orbital  function,  and  this  means  that  correlation  affects  the  electron 
distiribution  appreciably.  The  energy  drop  of  0.033  a.u.  «  0*90  eV  going 
fronn  the  bond  orbital  function  to  the  modified  electron-pair  fvinction  may  be 
regarded  as  resulting  from  this  intra-bond  correlation. 

The  calculated  value  of  the  dipole  moment  is  close  to  that  obtained 
by  tbe  SC  LCAO  method,  and  agreement  with  the  experimental  value  is  not 
excellent  but  reasonable. 

4) 

Comparing  these  results  with  those  for  water  one  cetn  see  that  their 
general  features  are  almost  the  same.  Thus,  we  may  conclude  that  the  success 
of  the  modified  electron -pair  method  is  bona  fide*  The  change  of  the  basic 
atonxic  orbitals  -  from  Slater  orbitals  to  SCF  -  and  the  change  of  the  number 
of  bonds  and  lone-pairs  do  not  affect  the  results. 


ACKNOWLEDGEMENT 

The  bulk  of  the  calculations  reported  in  this  paper  was  carried  out  on 
the  F^C  1  computer  at  the  Faculty  of  Science,  University  of  Tokyo.  The  authors 
are  grateful  to  the  members  of  the  Takahashi  Laboratory  for  help  in  using  the 
macbine  • 

Their  thanks  are  due  to  Professor  Kaplan  for  his  personal  communica¬ 
tion  concerning  typographical  errors  in  the  value  of  integrals  given  in  his  paper. 
They'  also  thank  him  for  the  values  of  the  dipole  moment  integrals  which  they 
used  • 

One  of  the  authors  (K.O.)  would  like  to  thank  Professor  Ldwdin,  the 
Director  of  the  Uppsala  Quantum  Chemistry  Group,  for  his  generous  hospitality 
and  Sill  the  members  of  the  group  for  the  benefit  of  a  stimulating  atmosphere. 

The  authors  are  much  indebted  to  Professor  M.  Kotani  and 
Dr.  J.W.  Eastman  for  reading  and  critisizing  the  manuscript. 


-1,0- 


REFERENCES 

1)  A.C<  Hurley.  J.E.  Lennard-Jones  end  J.A.  Pople:  Proc.  Roy.  Soc. 

A220  (1953).  446. 

2)  J.C.  Slateri  Quart.  Progr.  Rep.  Solid-State  and  Mol.  Theory  Group. 
Mass.  Inst.  Tech.  (15  April  1955).  unpublished. 

3)  J.M.  Parks  andR.G.  Parr  :  J.  Chem.  Phys.  2£  (1958).  335. 

4)  R.  McWeeny  and  K.  Ohno:  Proc.  Roy.  Soc.  A255  (I960),  367. 

5)  C.C.J.  Roothaaui:  Revs.  Modern  Phys.  ^  (1951),  69. 

6)  G.G.  Hall:  Proc.  Roy.  Soc.  A205  (1951),  541. 

7)  H.  Kaplan:  J.  Chem.  Phys.  26  (1957).  1704. 

8)  D.P.  Merrifield:  Quart.  Progr.  Rep.  Solid-State  ud  Mol.  Theory  Group, 
Mass.  Inst.  Tech.  (15  January,  1962),  unpublished. 

9)  P.O.  Lbwdin:  J.  Chem.  Phys.  ^  (1950),  365. 

10)  B.C.  Carlson  and  J.M.  Keller:  Phys.  Rev.  105  (1957),  102. 

11)  H.  Kaplan:  Private  commimication  to  the  authors. 

12)  G.G.  Hall:  Rep.  Progr.  Phys.  ^(1959),  1. 


z 


Figure  1 


The  co-ordinate  eyetem. 
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